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Genetic Algorithm for Lattice Gauge Theory. 

On SU(2) and U(l) on 4 dimensional lattice, how to hitchhike to thermal equilibrium state . 
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Applying Genetic Algorithm for the Lattice Gauge Theory is formed to be an effective method to minimize 
the action of gauge field on a lattice. In 4 dimensions, the critical point and the Wilson loop behaviour of SU(2) 
• lattice gauge theory as well as the phase transition of U(l) theory have been studied. The proper coding methodi 
| has been developed in order to avoid the increase of necessary memory and the overload of calculation for Genetic 
^\ ■ Algorithm. How hichhikers toward equlibrium appear against kidnappers is clarified. 
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The goal of numerical simulation is to get the 
maximum (minimum) value of a given objective 
function parametrized by some set of indepen- 
dent variables. A genetic algorithm originated by 
J.H.Holland Q is known that it is good at search- 
ing the maximum/minimum value under compli- 
cated situations. It bases its search on a set of 
point /individual of a searching space and on se- 
lection with respect to its fitness value, that is, 
the search of genetic algorithm is a global one 
rather than those of Metropolis and Heatbath 
method. Some severe extra works for coding a 
problem toward an appropriate form and treating 
some number of individuals at once had made its 
application wait until operating ability of com- 
puter system becomes adequate. By the recent 
development of the computer systems, its wide 
implementation becomes possible. However, it is 
still a subtle point that how to code a genome and 
how to design operations for the sake of good per- 
formance . We develop the coding of genomes in 
the SU(N) and U(N) lattice gauge theories and 
investigate to tune its searching parameters. 

1. Algorithm and Coding 

A configuration on a lattice is treated as a phe- 
notype. Link variables on a site are taken from 
the link variable pool created in advance by the 
integer number index which is coded to a gene, a 
binary string. In this coding, the memory for ge- 
netic algorithm operation is able to be suppressed 
as follows; With the number of link variables in 
a link pool, 2 l p, the lattice size, 2 Lx+L y +Lz+Lt , 
and the population size, 2 pop - stze 1 the length of a 
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genome becomes Lp x n_dimx lattice size. 

For instance, for 8 SU(2) lattice, a genome( a 
lattice) needs 282kb, and a population (32 lat- 
tices) needs about 9.2Mb. On the other hand, 
but for the genetic algorithm, the same size lattice 
needs 1.5 Mb and 32 lattices would need 48Mb. 
That is, a genetic operation needs only a quarter 
of memory for other methods. 

Our algorithm has three genetic operations; be- 
sides ordinal ones of selection and recombination, 
a new one we call "the intergenerational conflict" 
is introduced to establish the thermal equilib- 
rium. Furthermore as a local search, Metropolis 
updating applied to a link variable on a lattice 
is added, with which our combined algorithm be- 
haves like a hybrid genetic algorithm. 

This extraneous operation, the intergenera- 
tional conflict, is applied to the parent with a 
smaller action between a pair of parents' genomes 
and the offspring with a smaller action between 
the pair of offsprings produced from the parents. 
Offspring is passed to the next generation un- 
der the condition expressed using an ordinary 
Metropolis function, 

min{l,exp~ l3 ^ s ''~ Sp ^} > £, 

where 

S = S[U] = /3E P (\ - ^ReTrU^j , 

£ an uniform random number. 

Otherwise the parent is passed. Metropolis up- 
dating is applied to the survivor of the intergen- 
erational conflict. 

All genetic operations are coded in bit calcu- 
lations. Arithmetic operation which needs more 
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32*32 lattice beta=8.0 population size 128 
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32*32 lattice beta=8.0 population size 128 
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The horizontal axis is computing time. 
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The horizontal axis is the number of generation. 



Figure 1. Thermalizaion of Action per Plaquette 32 • 32 lattice (3 — 8.0 , population size 128 
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Figure 2. The variance of fitness, 32 ■ 32 lattice 
@ = 8.0 , population size 128 

executed time is used for the calculation of the 
action. This coding can reduce the increasing of 
the executing time. For instance, genetic opera- 
tions for 8 4 lattice with (3 = 2 occupies 17 per- 
cent of the executing time for one cycle, that is 
the increasing time caused by genetic operations 
is 15 percent of Metropolis method. This analysis 
conditions our genetic algorithm to an advantage 
against the usual Metropolis method. 



2. The algorithm dynamics 

To discuss the profitable condition, the dynam- 
ics of our algorithm on the statistical point of view 
is analized. For the procedure of Metropolis up- 
dating on a lattice is a random selection of a link 
=>■ decoding information of gene =>■ fixing the in- 
dex of the link variable =>• updating the variable. 
This Metropolis updating makes configurations 
converge fast or slow and or not converge at all, 
because it work not only on the target lattice but 
also on the other lattice in the stack. This up- 
dating duplication makes two cases occur. First 
one is that if Metropolis updating works well even 
for another lattice configuration, the other lat- 
tice's action decreases. That is, its fitness value 
increases, and the lattice which is not the target 
of the Metropolis updating does hitchhike. Sec- 
ond one is the opposite case, that is, if Metropolis 
updating works well for the target lattice config- 
uration but for the other lattice in the stack, its 
action increases that is its fitness value decreases. 
In this case, the target lattice kidnaps the other 
one. For the kidnapped lattice, the probability to 
be selected as a parent becomes low. 

The behaviour of Hitchhiker and Kidnaper are 
clarified from the statistical value of the fitness. 
Five different Metropolis updating ratios, from 
0.1 percents to 10 percents are set. From Fig. |l|, 
the Metropolis updating ratio is found to de- 
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Table 1 

U(l) Action on 4 dimension 
(3 Metropolis GA 
1.0 0.538 0.538 
1.5 0.370 0.369 
2.0 0.298 0.292 
3.0 0.237 0.235 

Table 2 

SU(2) Action on 4 dimension 
(3 Metropolis GA 
2.0 0.565 0.564 
4.0 0.384 0.34 
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Figure 3. Thermalization of the U(l) action on 
8 4 lattice j3 = 2.0, GA population size 16 and 
Metropolis method 

pend strongly on the performance, so that the 
Metropolis updating ratio should be set 1 percent 
and not 5 percent nor 0.5 percent. 

In Fig. |l| with generation, the effect of kid- 
napers and hitchhikers are appeared. There is not 
so clear difference between 10 percent line and 5 
percent line. The slowing down of 10 percent line 
is caused by the kidnapers. On the other hand, 
the line of 1 percent and 5 percent fall rather 
fast, which shows there are some hitchhikers in 
the populations. From the variance of fitness, the 
populations which have some hitchhikers keep the 
variance at some level not so high and not so low. 
It means that some lattices might be kidnapped 
but other hitchhiking lattices accelerate its ther- 
malization speed. It means the ratio of Metropo- 
lis updating has to be tuned. We set 0.01 for 
13 = 8.0 and 0.005 for (3 = 2.0 



3. Result and Discussion 

The physical values obtained by our algorithm 
are consistent with those given by Metropolis. 
See Tab. | and |. 

The mechanism Hitchhiker and Kidnaper is 
clarified from the statistical value of the fit- 
ness. Since the Metropolis updating ratio de- 
pends strongly on the performance, it should be 
set 1 percent not 5 percent nor 0.5 percent to 
(3 = 2.0. 

Without any tuning parameters, a performance 
of our algorithm for U(l) on 4 dimensional lattice 
with 16 population size and [3 = 2.0, is same as 
that of Metropolis method. See Fig. [| 

For SU(2), on 2 dimensional lattice, our algo- 
rithm can get higher performance than Metropo- 
lis method 0]. On 4 dimensional lattice, how- 
ever, it's performance is about the same as that 
of Metropolis. We have not tuned the parame- 
ter for hitchhikers on 4 dimensional lattice yet, 
but the performance is strongly dependent on 
the Metropolis updating ratio, because it occu- 
pies calculation time mainly. 

Our work suggests that an algorithm with ge- 
netic operation does work in the case that the 
thermal equilibrium state is to be established. 

The introduction of fermion must follow as the 
next stage, since our coding is good for SU(3) 
lattice gauge theory. 
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